Loose mechanochemical coupling of molecular motors 
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In living cells, molecular motors convert chemical energy into mechanical work. Its thermodynamic 
energy efficiency, i.e. the ratio of output mechanical work to input chemical energy, is usually high. 
However, using two-state models, we found the motion of molecular motors is loosely coupled to 
the chemical cycle. Only part of the input energy can be converted into mechanical work. Others is 
dissipated into environment during substeps without contributions to the macro scale unidirectional 
movement. 
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Introduction. In biological cells, molecular motors 
are individual protein molecules that are responsible for 
many of the biophysical functions of the cellular move- 
ment and mechanics They work in nanometer range 
and convert chemical energy, stored in ATP molecules, 
into mechanical work. Important examples of molecular 
motors are kinesin, dynein Q and mysion 0, Recent 
experimental data indicate that the energy efficiency of 
molecular motors is high. 

In literatures, there are many theoretical models to 
study molecular motors, such as Fokker-Planck equa- 
tion Langevin equation lattice model Q, net- 
work model [8j and etc. All of the existing models can 
be roughly classified into two categories: (1) There is 
only one chemical state, or equivalently only one tilted 
periodic potential in the model. (2) There are multiple 
chemical states, i.e. there are several periodic potentials 
in the model. To some extents, any one model of the sec- 
ond class is equivalent to one of the first class [9] . Usually 
the first class models are simple, and can be used to get 
explicit formulations of important biophysical quantities, 
such as the mean velocity [TO, [H| , effective diffusion co- 
efficient [l2| and mean first passage time [l3j . However, 
due to the oversimplification, more reasonable biophysi- 
cal properties of molecular motors can not be obtained 
from them. On the other hand, though the second class 
models seem more reasonable, it is too difficult to obtain 
explicit results [3, [H| . The simplest model of class two 
is the two-state model, it has been frequently employed 
by many authors jl6l4l9j. But most of their results are 
based on numerical calculations and then difficult to do 
further analysis. In this research, we will give some ex- 
plicit results of the two-state model, including the mean 
velocity and energy efficiency. The main conclusion that 
we draw from the two-state model is that, the motion of 
molecular motors is loosely coupled to the chemical cy- 
cle, its mean velocity might be zero even if there exists 
nonzero input energy, which is in accordance with the 
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where f(t) is Gaussian white noise, £ is viscous friction 
coefficient, fc# is Boltzmann's constant, T is the absolute 
temperature, and $(x) is a tilted periodic potential, i.e. 
$>(x) — $>(x + L) = A<£ is constant (the step size L = 8 for 
conventional kinesin and cytoplasmic dynein, and L = 36 
for Myosin V). Or equivalently, it can be modeled by the 
following Fokker-Planck equation 



dtp = d x (d x ^p/^ + Dd x p) . 
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where D = ksT/^ is free diffusion coefficient, p{x,t) is 
the probability density of finding the motor at time t and 
position x. Using either (pp) or (|2j), the steady state mean 
velocity can be obtained as follows 
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with f3 = l/k B T. Obviously V > if and only if A$ > 
0. If there exists external load F extl the potential $(x) 
should be replaced with <j)(x) = $(x) — F ext x. 

Another popular one-state model is the one- 
dimensional hopping model: in which the motor in 
mechanochemical state j can jump forward to state j + 1 
with rate Uj, or backward to state j — 1 with rate Wj. 
After moving N states forward or backward the motor 
comes to the same chemical state but spatially shifted 
by a step size L. Using this one-dimensional hopping 
model, the mean velocity of the molecular motors is 
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Meanwhile, the free energy difference in one 
mechanochemical period is Ap = ksTln (J^JSq 1 wi/u^j. 
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FIG. 1: Schematic depiction of the two- state lattice models: 
(Up) general case, (Down) special case, in which uf a — u % d — 
for i ^ M, N. 



Obviously Vn > if and only if A/i > 0. Actually, it can 
be proved mathematically that liniAr^oo Vn = V and 
A/i = A<3> (A/i = if there exists external force, see 
|24j). By (j3j) or (j4|, the stall force of molecular motor can 
be easily obtained F s = A<£/L, and the thermodynamic 
energy efficiency is rj = F ext L/ A<£, which increases with 
the external load F extl and r](F s ) = 1. Therefore, the 
one-state models are tightly mechanochemical coupled 
models. The maximum of energy efficiency is 1, which is 
attained at stall force F s . 

Two-state models. The general two-state contin- 
uous model is the following mechanochemical coupled 
Fokker-Planck equations: 



d t P =Dd x (f3<f>' 1 P + d x P) 
d tP =Dd x (P& 2 p + d xP )- 



hoj d (x)p-uj a (x)P 
Ud(x)p + u a (x)P, 



(5) 



with < x < L. Where P(x, t) and p(x, t) are the proba- 
bility densities of finding motor at time t, position x and 
in states 1 and 2 respectively. uo d (x)^uj a {x) are transi- 
tion rates between the two chemical states at position x. 
The corresponding general two-state lattice model is as 
follows (see Fig. [I]) 
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with n = 1,2, ••• ,iV, P n and p n are probabili- 
ties in states 1 and 2 respectively. All the rates, 
uj^uj^ are periodic with period N. In the 
following, we will only discuss the two-state lattice model 
(j6j), but all the corresponding results also can be obtained 
by the two-state continuous model (J5j) [25j. 

Although the steady state solutions of the general two- 
state lattice model (j6j) can be obtained explicitly, it will 
be certainly easier for us to show our conclusions only 
by employing one special case, in which uj l a — uj l d = 



for i ^ M, N. For the simplicity of notations, we denote 
0J^,uJ d by ft a ^d and uj^ 4 \ujf by u a ,u) d respectively 
(see Fig. [I]). To this special case, one can calculate that 
the probability flux is 
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Similar expressions for r k , s k ,t k , g kl h k can be ob- 
tained by replacing Fi^Bi^Vta^Vtd in expressions of 
R k ,S k ,T k ,G k ,H k with fi,bi,Q d ,Q a respectively. Under 
no external load, the input energy in unit time is 
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One can show that A/i = if and only if Rn = rN = 1 
and J c = 0. Rn = = 1 means the potentials in states 
1 and 2 are all periodic and so no input energy in either 
of the two states. J c = means there is no input energy 
during transitions between the two states. 

In the following, we assume that, under no external 
load the potentials in each states are all periodic since, 
biophysically the input energy to molecular motors comes 
from ATP molecules, which is hydrolyzed during transi- 
tions between different chemical states. So, under no 
external load, the mean velocity of the motor is 

V = JL = (Sn/Tn - s N /t N ) JcLPn/W, (9) 

and the input energy in unit time is 

A/i = k B TJ c P N [\n(oj a n d RM/tt a uJdrM)} /W. (10) 

Obviously, A/i = implies V = 0, but V = does 
not give A/i = 0. The reason is that there might exist 
nonzero circumfluence of probability (see Fig. [3]), i.e. the 
motors might make a forward substep in one state but 
stepping back in another one. In these processes, free 
energy is consumed but no effective steps are made. 

under nonzero external load F ex t, the tran- 
sition rates satisfy Fi(F ext ) / B i+1 (F ext ) = 
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[Fi(0)/B i+1 (0)]&cp(-5lpF ext L), fi(F ext )/b i+1 (F ext ) = 

[/i(0)/6i + i(0)]exp(-5^F esct L), in which ^Jo' 4 = 1 
(j = 1,2). Depending on the parameters Q a , cj^, cj a , f^, 
the motor might move to the right direction (V > 0) 
or to the left direction (V < 0). In our discussion, we 
always assume that the direction of the external load is 
opposite to the motor motion direction, i.e. the external 
load points to the left (F ext > 0) if V > and points to 
the right (F ext < 0) if V < 0. This assumption means 
F ext V > 0. 

To the nonzero external load cases, the free energy 
difference in one mechanical cycle is [25| 

AG =k B T{[(R N - 1)W/T N - J C (T N - S N )/T N ]]nR N 
+ [(r N - l)U/t N + J c (t N - s N )/t N ] In tat 
+ J c ln[(u a n d R M )/(ttaUJdrM)}} Pn/V. 

(11) 

It can be verified that the input energy A[i(F ext ) = 
AG(F ext ) + F ext V. In which AG(F ext ) > 0, and 
AG(F ext ) — if and only if J c = 0. Therefore, 
AG(F ext ) > if and only if Afi(F ext ) = 0, which means 
FextV < Afi(F ext ) and the equality holds if and only if 
FextV = Afi(F ext ) = 0. 

The velocity-force relation can be obtained by formu- 
lation ([7]) with the stall force F s is attained when J = 0. 
It should be pointed out that, under stall force, V = 
but the input energy Aji(F ext ) ^ if J c ^ 0, which 
only depends on the parameters fi a , uod-, oj a , Qd- In fact, 
J c = [u a Q d RM(0) - Q a Udr M (0)}exp(-Sl3F ext L) with 
S = S J + • • • + S J M _ 1 (j = 1,2). Therefore, the motion of 
molecular motors is loosely coupled to the chemical cy- 
cle, this is consistent with the recent experimental results 

[lill. 

The flashing rachet model can be regarded as one of 
the special cases of the above models: in which one of 
the potentials (for example, in state 2) is constant, i.e. 
/i(0) = &i(0) = / (i = 1, 2, • • • , N). To this much special 
case, corresponding results can be derived easily. 

Limit properties and energy efficiency. To 
some molecular motors, it might be possible to 
change the transition rates Q a , uod, uo ai f^, which 
usually depend on temperature, ATP concentration 
and potential profiles. For simplicity, we suppose 
(O a , uod, cj a , ^)=A(1^, (j®, Qj) and A is a parameter 
that can be changed experimentally. 

One can easily show that the mean velocity V — » 
with A — » 0. But for A — >• oo, the mean velocity tends to 
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Similarly, as A — > oo, the input 

energy Aji tends to /c#T\I/i \j\.{k\Rm / ^2^m) So, the 
output power of molecular motor has a limit which only 
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FIG. 2: The thermodynamic energy efficiency of molec- 
ular motor calculated by the special two-state model, 
as a function of external load F ex t (left) and as a 
function of parameter A (right), which is defined by 
(O a , cjd, cj a , Hd)=A(Qa? 0;°, &>d)- The parameters used in 
the model are N = 2, M = l,k B T 1,F (0) 100, Si (0) = 
0.5,Fi(0) = 0.005, B o (0) = 1,/ (0) = 0.5,&i(0) = l,/i(0) = 
50, b (0) = 25, L = 1, Q a = 20, Q d = 5, u a = 10, u; d = 50, 8{ = 



0.9,^=0.1 0=1,2). 



FIG. 3: Different cases of probability flux in which there are 
only two nonzero transition rates. 



depends on the parameters ki,k,2 (i-e., is independent of 
A). 

Another interesting biophysical property of molecular 
motor is its energy efficiency. By ([7]) (fTO]) . its thermody- 
namic energy efficiency r](F ext ) := FV/Afi is 
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Obviously, 77(0) = rj(F s ) = 0. From the loose 
mechanochemical coupling discussion in the above sec- 
tion, we can easily know 77 < 1. With F ext increases from 
zero to F s , the mean velocity decreases to zero monoton- 
ically. But for the efficiency 77, there exists a maximum 
value between and F s (see Fig. [2}. 

Discussions. The analysis of two-state models indi- 
cates the motor motion is loosely coupled to the chemical 
cycle. The biophysical reason of this loose mechanochem- 
ical coupling is there exist forward and backward sub- 
steps, which had been found experimentally in [26|, [27j. 
Free energy is consumed in each substeps, but only part 
of them really contributes to the effective unidirectional 
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FIG. 4: Two-state model of motor protein kinesin: In state 1, 
the motor moves from location A or A' to the potential well 
B, in state 2, the motor moves stochastically from location 
B to A or A' . The mechanical stepsize [A'A\ is about 8 nm, 
and the substep \A' B\ is about 5 nm }28l ]. 

motion. In the two-state models, this corresponds to the 
existence of probability circumfluence. Due to the tem- 
poral symmetry of states 1 and 2, there are altogether 9 
different cases of the probability flux (see Fig. [3]). Free 
energy is consumed but without effective forward motion 
in each circumfluence. To the general two- state lattice 
models, there exist much more cases of the probability 
circumfluence. 

Using special cases of the two-state continuous 
models ((5]), the same conclusion can be obtained. For 



example, to the special case: uo a {x) = 0Jd(x) = for 

< x < a,a < x < L, the explicit expression of 
the velocity gives: the mean velocity V > if and 
only if [n a cj d e^(^( )-^(°))-a; a n d e^^( )-^( a ))] x 

J L e^ 1 ^ // Q a e^ 1 ^ - j L e^ 2 ^ /J* e^ 2 ^' < 0. 

But the input energy is always positive. Therefore, the 
fact of loose mechanochemical coupling also can be found 
from this continuous models. As an example, the motion 
of motor protein kinesin can be schematically described 
by the special case of our two-state model (see Fig. H]). 
In which the transition from state 1 to state 2 is due 
to the neck linker docking of the microtubule bounded 
head (leading head), the transition from state 2 to state 

1 is due to ATP hydrolysis and phosphate release. Due 
to the existence of substeps, which might have no any 
contributions to the macroscopic unidirectional motion 
of kinesin, the energy efficiency of kinesin is certainly 
smaller than 1. 
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